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Abstract 



We study renormalization group flows in the Lifshitz-like iV-flavour four 
fermi model discussed in 0905.2928. In the large- N limit, a nontrivial flow 
occurs in only one of all possible marginal couplings and one relevant cou- 
pling, which provides the scale for Lorentz invariance violations. We discuss 
in detail the phase diagram and RG flows in the space of couplings, which 
includes the Lifshitz fixed point, the free field fixed point and a new fixed 
point characterized by z = 1 scaling and a violation of Lorentz invariance, 
which cannot be tuned away by adjusting a parameter. In the broken sym- 
metry phase, the model flows from the z = 3 Lifshitz-like fixed point in the 
ultraviolet to this new fixed point in the infrared. However, in a modified 
version of the present model, which has an effective ultraviolet cut-off much 
smaller than the Lorentz invariance violating scale, the infrared behaviour is 
governed by an approximately Lorentz invariant theory, similar to the low 
energy limit of the usual relativistic Nambu— Jona-Lasinio model. Such a 
modified model could be realized by a supersymmetric version of the present 
model. 
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1 Introduction 

Recently, there has been considerable interest in field theories in which 
Lorentz invariance is explicitly violated by terms containing higher order 
spatial derivatives. The presence of these terms leads to softer ultraviolet 
behaviour while preserving unitarity (which is typically lost in the presence 
of Lorentz invariant higher derivative terms because of ghosts associated 
with higher time derivatives). At very high energies, these Lorentz invari- 
ance violating terms dominate, leading to Lifshitz-like anisotropic scaling 
symmetry (in the classical theory) in which time and space scale differently: 
x — > x/a,t — > t/a z . The exponent z characterizes the scaling symmetry. 

Lifshitz-like field theories with anisotropic scaling have been used in con- 
densed matter systems to describe quantum criticality [l]-[6]. Recently they 
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have also been discussed in string theory in the context of possible appli- 
cations of AdS/CFT duality [I]-[T5] to condensed matter systems involving 
strongly interacting constituents. In a separate development, the idea that 
a relativistic theory at low energies may have a Lorentz non-invariant ultra- 
violet completion was suggested in [TO]. This idea has been further explored 
in [UJ-[25]. The suggestion that an ultraviolet completion of quantum grav- 
ity may be similarly formulated [26] has serious difficulties [27]- [33] because 
gravity elevates Lorentz invariance to a local gauge symmetry, which cannot 
be broken except by some kind of Higgs mechanism. In the present work we 
will focus only on non-gravitational theories. 

Lifshitz-like field theories with Lorentz invariance violations (LIV) have 
also recently been discussed in the context of applications to particle physics 
|33l E51 HSl ESI EZ]. In [33 OS] it was argued that a z = 3 Lifshitz- 
like ultraviolet completion of the Nambu— Jona-Lasinio (NJL) model [38] 
in 3 + 1 space-time dimensions has the required properties to replace the 
Higgs sector of the electro-weak theory. The four-fermion coupling in this 
model is asymptotically free, leading to dynamical mass generation for the 
fermions and chiral symmetry breaking. In an appropriately gauged version, 
fluctuations of the magnitude of the fermion bilinear order parametei0 can be 
interpreted as the Higgs field. This model obviates the need for fine tuning 
at the expense of introducing LIV at high energies. The hope is that for a 
sufficiently large LIV scale, the low-energy theory could be consistent with 
experimental constraints on This may, however, require a new fine 

tuning of parameters [2T] . 

The main purpose of this work is to analyze renormalization group (RG) 
flows in the model of [19] to understand in detail the possible emergence of 
Lorentz invariance at low energies. The analysis has been performed in the 
leading large- N approximation. Our findings can be summarized as follows. 

• At low energies, in the fermionic sector, the theory recovers approximate 
Lorentz invariance, violations being of order E 2 //j, 2 , where [i is the energy 
scale associated with LIV. However, in the bosonic sector, in the broken 
chiral symmetry phase, the induced kinetic terms violate Lorentz invariance 

lr The phase of the fermion bilinear is the Goldstone mode which combines with the 
gauge field as usual to make it massive. 

2 For current situation on experimental searches of Lorentz symmetry violations, see 

[H-iii. 
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at 0(1) level. The origin of these violations is simple to understand — they 
arise from fermion modes with energies higher than \x propagating in a loop. 
These violations can be made small by imposing an effective cut-off on the 
theory and arranging fi to be much larger than the cut-off (see the last 
paragraph of Section [3]) • 

• The RG flows reveal a new nontrivial fixed point, apart from the z = 3 
ultraviolet fixed point. The theory flows down from high energies to this fixed 
point at low energies. The above mentioned approximate Lorentz invariance 
in the fermionic sector and 0(1) LIV in the bosonic sector are characteristic 
of this new fixed point. If one works with a fixed finite cut-off A, and sends the 
LIV scale fi much above A, the Lorentz violations become smaller and smaller, 
leading to an approximately Lorentz invariant theory which is identical to an 
effective theory derived from the NJL model at low energies, with a cut-off 
A. 

The plan of the paper is as follows. In section 2 we deform the z=3 
action with all possible marginal and relevant couplings (from the point of 
view of z = 3 scaling) and study their effect on the vacuum solution. We find 
that only one of the three possible four-fermion (marginal) couplings has a 
nontrivial flow. Moreover, there is only one relevant coupling that affects the 
low energy physics, namely the coupling that determines the scale of LIV. 
This had already been remarked in [19], but section 2 provides a detailed 
justification for it. In section 3, we derive the effective action at a scale much 
smaller than the scale of LIV. We find that at low energies, in the broken 
symmetry phase, the kinetic terms induced by fermion loops for the massive 
scalar bound state of the fermions violate Lorentz invariance at 0(1) level, 
which cannot be corrected by any fine-tuning of parameters. In section 4 we 
study the RG flow of the couplings and locate fixed points. We compare this 
fixed point structure with that in the relativistic NJL model. We end with 
some concluding remarks in section 5. Details of some of the calculations 
have been given in three Appendices. 

2 The action and vacuum solutions 

The model discussed in [19] consists of 2N species of fermions, ip a i(t, x), where 
the index a runs over the values 1, 2 and the index i runs from 1 to N. Each 
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of these fermions is an SU(2) spinor, where SU (2) is the double cover of the 
spatial rotation group 50(3). It is useful to view the index a as denoting 
the two Weyl components of a Dirac fermion in a four dimensional theory 
with Lorentz invariance. The action we consider will have the following 
symmetries: a global U(l)i x Z7 (1)2 Symmetrel under which the fermions 
transform as 

i>ai ^ e ia «i> ai , a = 1,2; (1) 

and a global U(N) symmetry under which the fermions transform in the 
fundamental representation 

i>ai -> Uijlpaj, i = 1, N. (2) 

In addition to these symmetries, we will ensure that the action is invariant 
under the interchange ^u(t,x) — > ij}2i(t,—x). This is the analogue of the 
parity operation in the relativistic Dirac theory. 

A general action which is consistent with the above symmetries and con- 
tains all the relevant and marginal couplings is given byj 

S=Jd 3 x dt [id t - id.a (g Q (-id) 2 + g^j - g 2 d 2 } ^ 

+^2i [id t + id. a [g {-idf + - g 2 3 2 } i> 2i + g 3 {^u^u + ^ 2i ip2i) 



+9l \ {^li^li) 2 + (^2>2i) 2 \ + 9l tyli&li^ijfoj) + 9 2 (?P\i^2i^ 2j lplj) 



(3) 



At high energies, this action has the scale invariance of a z = 3 Lifshitz-like 
theory in which space and time have the dimensions [x] = — 1 and [t] = —3. 
At low energies, the relevant term with coupling g\ dominates, so one expects 
the model to flow down to an approximately Lorentz invariant theory at low 
energies. We will study the dynamics of this action in the large N limit. 



3 This is the analogue of chiral symmetry in the corresponding Lorentz invariant four 
fermi model. 

4 A11 other possible four-fermion terms, like e.g. {ip[i^2j){i>2 3 iV'ii): can be related to the 
terms in (|3|) and/or terms involving vector bilinears like (ip^ai/iii) .(ipljijip2j) ■ The latter 
do not affect vacuum solutions since the vcvs of such bilinears vanish because of invariance 
under spatial rotations. 
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It is useful to employ the notation of Dirac matrices and rewrite the action 
in terms of the four-component spinor 



(4) 



and its Weyl components P L = ^ = ijji,2i- We will also use the Dirac 
gamma matrices 

7° = a x <g> 1, Y = ia 2 ® a u 7 5 = iWtV = a 3 ®l, (5) 

where the a's are standard Pauli matrices satisfying = ieijk&k and 1 

is the 2x2 identity matrix. In terms of these matrices, the Weyl projection 
operators are given by Pl,r — |(1 ± 7 5 )- In this notation, the above action 
takes the form 



,5' 



d x dt 



{ 7 ido ~ 9id l +g 3 )+ il -d gi - g d 



:T2 



-g 2 {^Li^Ro){^R^Li 



(6) 



As is usual for actions with four-fermion interactions, we now introduce 
auxiliary scalar fields to rewrite the above in the completely equivalent form 
of an action quadratic in fermions: 



,5' 



d 3 x dt 



-7° (aP L + pP R ) Ui 



*A 7° (ido ~ g 2 d 2 +g 3 )+ ird U - god 2 ) + {<\>P L + <P*Pr 



p 2 + pv gffi 



g\ 



gi 



g* 



(7) 



Here p and rj are real scalar fields and is a complex scalar field. We have 
also defined 



2p + t) = a, 



9 4p + 2 9 4v = P. 

gi 



(8) 



The scalars are fermion bilinear composites, as can be easily derived from 
their equations of motion: 



P = gl(^Li7°^u 



r) = gl^ml^Ri, 



g 2 {^ L ^Ri). 



(9) 
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2.1 Vacuum solutions 



Since the action ([7j) is quadratic in fermions, one can integrate these out 
to get an effective action for the scalar fields. For vacuum solutions, it is 
sufficient to consider only the homogeneous modes, po, rjo and 0o, of the 
scalar fields. Moreover, without loss of generality, one may take to be 
real. In this case, the effective action for the scalars is given by 



where V denotes the volume of space-time. Also, in the momentum space 
representation, we have 

l (k ,k) = k + g 2 k 2 + g 3 , l(k) = k{g k 2 + 9l ) , k — \k\. (11) 

The equations of motion are obtained by varying this action with respect 
to Poj Vo an d 0o- The calculation can be simplified by noticing that due to 
the "parity" symmetry ^ii{t, x) — > ^Ri(t, —x), which we assume is unbroken, 
and the relations (Q, the vacuum solutions po an d ?7o are related, i.e. 



Thus, there are really only two independent variables, po and 4>q . Using this 
relation after varying (TTUj) with respect to p , rjo and O > the two independent 
equations of motion can be written as 




-i N Tr In j/ o7 - l.j + O + tVo^l + PqPr) 
( P 2 o + PoVo , glvo . <Po\ T/ 



(10) 



po/gl = vo/gl 



(12) 



Po 
Ai 




d 4 k (Z + «o) 

(2vr) 4 (/ o + « o )2_r2_ 2' 



(13) 




(14) 



where we have defined the 't Hooft couplings, 



#4,5^ = A 4 ,5, 



g 2 N = A, 



(15) 



and made use of the relation 




(16) 
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which follows from equations flS]) and (TT21) . The equations (TT5|> and <HM are 
exact in the usual large N limit in which the 't Hooft couplings are held fixed. 

Let us first consider equation ffl~3l) . The momentum integral on the right 
hand side has a power divergence, so the dependence on a$ cannot be removed 
by shifting k . To do the integral, we need to introduce a regulator. We will 
use a simple cut-off regulator Continuing to Euclidean momenta, we may 
then write 



where I = \l\ = k\g k 2 + gi\, g{k) = (g 2 k 2 + g 3 + do) and p = ipo, 92 = igi 
and g 3 = ig 3 are Euclidean continuations of respectively p , g 2 and g 3 . 

In the limit of large M, the integral in (ITT)) can be done by expanding 
both the numerator and denominator of the argument of the 'In' function 
around (M 6 + g^k 5 ). Discarding terms that vanish as M — > oo, we get 



Po = ~i \ ~M 2 g 2 I - (g 3 + a )h + (2g o9l - gf) g 2 I 2 + -glh \, (18) 



where the coefficients Iq, Ji, I 2 and ^3 are functions of go only and are 
listed in Appendix 1. In writing the above, we have continued the various 
parameters and fields back to Minkowski signature. Now, using relation f|T6|) . 
the above equation can be solved for p : 



Notice that po vanishes if g 2 = g 3 = 0. If g 2 — but g 3 7^ 0, then 
Po is nonzero, but we may take all the couplings appearing in the solution 

5 Note that the regulator used here explicitly breaks Lorentz invariance. This does not 
affect vacuum solutions. However, in sections 3 and 4, where we consider the question of 
restoration of Lorentz invariance at low energies, it will be important to choose a suitably 
different regulator that allows for such a possibility. 




(17) 
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as fixed, independent of the cut-off. In case g 2 7^ 0, then the right hand 
side of (1191) diverges quadratically with the cut-off. However, this divergence 
can be removed by shifting g 3 , i.e. we set g 3 = (g' 3 — M 2 g 2 Io/ where 
g' 3 is independent of the cut-off M. This determines po i n terms of the M- 
independent couplings g 2 , g' 3 , g$, etc. 

The other vacuum parameter, O , is determined by equation (fl4l) . In 
this case, the integral over ko on the right hand side is convergent and the 
entire integral is only logarithmically divergent. So a shift of the integration 
variable k is allowed. Doing this enables us to get rid of the couplings 
92, 93, 9a, 95 and po from this gap equation. After a Euclidean continuation 
(ko = ik A ) it then takes the form 

\_ !? [d^k 1 

x~ J ^(kjTmTWr ( } 

As discussed in [19], this equation signals the breaking of the global U(l)i x 
f/(l)2 "crural" symmetry (OQ) down to 17(1). It also requires that the coupling 
A have a nontrivial RG flow, which depends on the couplings go and g\. We 
will analyze this flow in detail in section 4. 

We end this section with the following comment. It should be clear from 
the above discussion that in the present model all the interesting dynamics 
arises from the NJL type of four fermi interaction (conjugate to the 't Hooft 
coupling A) and its RG flow is decoupled from the dynamics of all the other 
marginal couplings. Therefore, for calculational simplicity, in the rest of the 
paper we will set the marginal couplings g^ and g^ to zero. We will also set 
92 — 93 — since these couplings also do not affect the RG flow of A. 



3 Low energy effective action 



Setting g 2 = g 3 = #4 = g$ = in (jHJ) as discussed above, and going over to 
the bosonic variable 4>(x) as in (J7J), we get the action 



,5' 



d 4 x 



^ l <i-y°d + i^.d(g 1 - g d 



( < f ) p L + < f ) *P R )^ i -!^\ 



(21) 



We are interested in the low energy effective action for this system in the 
phase in which the symmetry (OQ) is broken. To derive the low energy action, 
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/ , . B(x) 

we substitute in the above <f>(x) = (0o + ^^) e% ^ i where 0o is the vacuum 
solution and cr(x) and 6{x) are respectively the magnitude and phase of the 
fluctuation of 4>(x) around the vacuum. We get, 



,5' 



d x 



wArfdo + 17.9 [gi - g d 



a 



P L e 1 ^ + P R e \% 



—\ G ' 



2VN 
A 



The phase 9 is the Goldstone model- It can be "rotated" away from the 

i e 

Yukawa- type coupling by making the change of variables —> e 

and — *> e' 2 ^^^. Its derivatives will, however, now appear from the 

fermion kinetic terms. After this change of variables, the action takes the 

form 



5 



d x 



i^do + ij.d ( 01 - goB 2 ) + O + 



a 



A 



-a 



2VN 
A 



N 



(23) 



where 0(^89^) indicates terms involving derivatives of the Goldstone mode 
9. In the following we will ignore these terms and retain only the mode a. 
There is no basic difficulty in retaining 9 also, but this is not essential for 
describing the RG flow of A and LIV properties of the low energy action^. 

To proceed further, it is convenient to rewrite the action (123]) 
turn space. We have 



m momen- 



5 



d A k 



d A q 



(2vr) 4 J (2tt) 



(22) 



(7% - Tl(k) + O ) (2nf5\k -q) + -La(k - g)Wg) 



(2tt) 4 A 



2VN 
A 



<h*r(k = 0) 



VN 



6 If the U(l)i x U(l)2 symmetry is gauged (in such a manner that there are no anoma- 
lies), then this mode will be eaten up by the gauge fields by the usual Higgs mechanism. 

7 There are no terms in the low energy action which mix 9 and a at the quadratic level. 
Moreover, the induced kinetic terms for both the fields are finite in the limit of infinite 
UV cut-off. Only the coefficient of the a 2 term diverges and needs to be renormalized. 
Hence, for simplicity, in the following we have set to a constant. 



(24) 



9 



where V denotes the volume of space-time and l(k) is given by (11 II) . As 
discussed in the next section, although this action depends on the two pa- 
rameters go and g%, the physics described by it actually depend on only a 
particular combination of these, namely 

^=|0i/0o|*. (25) 

ii has the dimensions of energy. As discussed below, its significance is that 
it is the energy scale at which Lorentz violations become important. The 
action (124)) can, in fact, be explicitly written in terms of this combination. 
However, two different versions of the action are needed to cover the entire 
line of real values of /a, from zero to infinity. In the form (1241) . the action is 
valid for all values, but at the expense of having a redundant variable. We 
will discuss this issue further in the next section after we have obtained the 
low energy action. 

The action at low energies is obtained by integrating out high energy 
modes from this action. To carry out this procedure, we will need to reg- 
ularize divergences in loop momentum integrals, which we will do by using 
a cut-off that scales like energy. However, the way in which the cut-off is 
imposed needs to be chosen carefully, because an arbitrarily chosen cut-off 
procedure (in a theory with LIV) might preclude the possibility of an approx- 
imate restoration of Lorentz invariance at low energies. We choose to impose 
a cut-off on the expression {k\ + l 2 (k)), where k± = —iko is the Euclidean 
continuation of energy. To understand why it is natural to adopt this cut-off 
procedure, recall the definition of l(k) = \l(k)\ = k\g Q k 2 + g x \, given in (TlTj) . 
Now, for gi ^ 0, we can completely scale it out of this expression by the scal- 
ing k — > k/\gi\. This gives / — > k\efi~ 2 k 2 + 1|, where e = ±1 is the relative 
sign of go with respect to g\. For k << ll, (k 2 + l 2 {k)) can be well approxi- 
mated by the 5'0(4)-invariant form {k\ + k 2 ). Thus, the proposed cut-off is 
naturally consistent with the possibility of Lorentz symmetry restoration at 
low energies^]. 

In the following we will use the notation f[d A k]E->A to indicate that the 
integral over Euclidean momenta is restricted to the region E 2 < (k\ + 

8 By contrast, a cut-off such as the one used in (|17[) is intrinsically Lorentz non-invariant . 
With such a cut-off, possible emergence of Lorentz invariance at low energies will not be 
manifest. 
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l 2 (k)) < A 2 , i.e. 



[d 4 fc] B _ A ••• = 

J d 4 k 9 (A 2 - (kj + l 2 (k))) 9 ((kl + l 2 (k)) — E 2 ) • • • , (26) 

where the 9 is the usual step function. Note that both A and E have dimen- 
sions of energy. For the special case of E = 0, it will be convenient to use 
the more compact notation j[d 4 k]\. 

Let us now integrate out from (1241) all the modes between {k\+l 2 {k)) = A 2 
and (k\ + l 2 {k)) = E 2 , where A is the cut-off and E < A. The action with 
the lower energy cut-off E is given by 



- J Jpjf (C„ + C,kl - C,e) M*)| 2 - 2C 3 ^ a(k = 0) - 

+ classical part. (27) 

Details of the calculation leading to fT2T|) as well as the calculations of the 
coefficients C n (n = 0,1,2,3) have been given in Appendix B. We have 
retained terms only up to Oil) in large N and up to two derivatives of a(x). 
The "classical part" refers to the contribution which comes from the action 
evaluated on the classical solution. 

Notice that even though the scalar field started out as an auxiliary field, 
it has now developed kinetic terms. As can be seen from equations fl56|) - (l58|) . 
these terms imply a maximum attainable velocity at low energies for the a 
particle, which is given by 



4 = | = 4 % (28) 



where = \gi\ is the maximum attainable velocity of the fermions at low 
energies. The quantities C*i i2 have been defined in (1581) . In general c a ^ c^, 
leading to LIV at low energies. As a check on our calculations, it is not 
difficult to see that for g = 0, c a = \gi\. Thus, we recover properties of the 
usual relativistic NJL model in the appropriate limit. Note that for go ^ 0, 
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the model is renormalizable since the momentum integrals involved in the 
computation of are finite in the limit A — > oo. However, for g = 0, a 
finite A is essential since in this usual, nonrenormalizable, relativistic NJL 
case the momentum integrals diverge. 

What is the extent of the LIV at low energies? It follows from the above 
discussion that the magnitude of LIV depends on the deviation of the ratio 
c a /c^ from unity. From equations (156|) - (!58|) in Appendix B, we see that this 
quantity depends on the various energy scales only through the dimensionless 
ratioa^l 

fj, R = fi/E, p, = fi/A. (29) 

In Figure [1] we have plotted c a /c^ as a function of /i^ 1 for different values 
of fl. In these plots, we have chosen the relative sign of go and g\ to be 




Figure 1: The variation of the ratio c a /c^, or equivalently (C 2 /Ci) 1 ^ 2 , cf. 
(128]) . with /i^ 1 for different values of fl. This ratio measures Lorentz invari- 
ance violations. The dashed curve corresponds to \x — 0. The other two cases 
have fl = 1 and 5. We have used (fto/n = 10 -3 . 

positive, i.e. e = +1. The dashed line in Figure [1] corresponds to case with 
fl = 0, i.e. the case in which the cut-off A is infinite. We see that even at 
very large values of /i# there is violation of Lorentz invariance at 0(1) level. 
The origin of these LIVs becomes clear from the other two lines in Figure 
[TJ The full line corresponds to the case with fl = 1 and for the dotted line 

9 The reason for the subscript "R" will be clear in the next section where we discuss 
the renormalization of the various couplings. 
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fl = 5. We see that the LIVs are dramatically reduced in these cases, being 
smaller for larger value of /L Clearly, fermion modes with energy larger than 
\i propagating in the loop (Eq. (150)) ) make Lorentz violating contributions to 
C12 at 0(1) level. These modes can be removed from the loop by imposing a 
cut-off which is smaller than the LIV scale, i.e. for p > 1. For smaller cut-off, 
the effect should be better. This is precisely what we see in the calculations 
shown in Figure [D 

It should be emphasized that the 0(1) violations of Lorentz invariance at 
low energies that we have found in the bosonic sector of the present model 
cannot be tuned away by adjusting any parameter. This makes (an appro- 
priately gauged version of) the present model an unsuitable candidate for an 
alternative solution to the hierarchy problem. A way out could be provided 
by a supersymmetric extension of the model. If supersymmetry is broken at 
a scale M s much smaller than the Lorentz invariance violating scale /1, then a 
cancellation with bosonic partners would remove the Lorentz violating con- 
tributions from fermionic modes with energy larger than M s propagating in 
the loop. This essentially means that the role of the cut-off A would then 
be played by M s . The residual low energy LIV would then be controlled 
by M 2 /yU 2 , which can be made small by tuning the scale at which super- 
symmetry is broken. Existing constraints on LIV coming from bounds on 
the maximum attainable velocity of various particles (6c < 10" 24 ) [39]- 02] 
implies M s /fi < 1CT 12 . For /i ~ 10 17 GeV, this means that supersymmetry 
can be broken at a much higher scale, M s ~ 100 TeV, than in the currently 
popular scenarios^. The present scenario might become a serious possibility, 
should LHC see a composite Higgs but no supersymmetry. 



4 RG flows and fixed points 

Classically, at high enough energies (where all masses can be ignored), the 
action in fl23|) has the scaling symmetry under which x° —>■ x°/a 3 , x — > 
x/a, ^ — > a 3//2 \I/, cr — » a 3 a, where a is the parameter of scaling transforma- 
tion. This is the z = 3 Lifshitz-like fixed point. For go = 0, classically the free 
fermi theory has the scaling symmetry x° — > x°/a, x — > x/a, ^ — > a 3//2 \l/. 
This is the familiar Lorentz invariant case, which can be described in the 

10 For a recent review, see [43] . 
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above language as a z = 1 fixed point. At energies much smaller than the 
scale [a, Lorentz violations are small even for z = 3, and so classically one 
recovers an approximately Lorentz invariant theory at low energies. In the 
following, we will study RG flow in the quantum theory from the z = 3 fixed 
point in the ultraviolet to find out what theory it flows to in the infrared. 

4.1 Determination of renormalized parameters 

In order to study the flow from high to low energies, we need to find out how 
the various couplings get renormalized. The starting point in the determi- 
nation of the renormalization of the couplings in action (|24|) . in the leading 
large N approximation, is the low energy action (|27|) . To implement the 
Wilson RG procedure, we need to rescale the cut-off E in (1271) back to the 
original cut-off A. As discussed above equation our cut-off procedure 
imposes the restriction (fc| + l 2 {k)) < A 2 on Euclidean space momentum in- 
tegrals. Writing E = bA, we see that the cut-off E on the low energy action 
( 1271) can be rescaled to A by the scale transformations (change of variables) 
A; 4 — > bki, k — > ak, followed by the scalings of the couplings 

a 3 9o = bg' , ag x = bg[. (30) 

These give the renormalized couplings and RG flow in the free fermion theory. 
Notice that the scaling parameter, b, for the energy is, a priori, unrelated to 
the scaling parameter, a, for the momenta. The z = 3 fixed point behaviour 
corresponds to choosing a = b 1 ^ 3 , and then the couplings scale as g' = g and 
g[ = b~ 2 / 3 g\. Since in this case go is invariant under the RG flow, we can set 
it to unity by scaling k — > |go|~ 1//3 &, leaving only one independent coupling, 
namely g\. Choosing a = b instead, one gets g' Q = b 2 ^ 3 g and g[ = g\, 
which are the scalings appropriate for a z = 1 fixed point. In this case, g\ 
is invariant under the RG flow and so we can scale it away by k — > k/\g\\. 
Once again we are left with only one independent coupling. 

The two fixed point behaviours discussed above can be treated together 
by setting a = b x l z where z = 3 or 1 0. Then, using the cut-off A to define 

11 Away from the fixed points, in general the RG equations (|30|) describe a flow in two 
parameters, namely a and b. More generally, in a theory with anisotropy in n different 
directions, the RG equations will describe an n-parameter flow, ft would be interesting to 
explore such more general flows. Here we will confine ourselves to a more traditional view 
of RG as a flow in a single scale parameter. 
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the dimensionless renormalized couplings, we get 

g 0R = Af" 1 ^ = E'-igo, g 1R = A^Vi = E^ 1 g 1 . (31) 



They satisfy the RG equations 



9or = - ( 1 J 9or, 9ir= - ( - - 1 ) 9m, (32) 



z 



where a dot denotes a derivative with respect to (— \nE) 0. Using these two 
couplings, we can define the renormalized version of (j25|) : 

\9ir/9or\* = t*/E = t*R- ( 33 ) 

This is precisely the quantity defined in (|29|) . 

In the leading large- N approximation, the free field renormalization (J32l) 
is not affected by the Yukawa coupling. However, the 't Hooft coupling A 
does receive quantum corrections. Its renormalization can be deduced from 
the term proportional to Co in the low energy action (]27p . Scaling the cut-off 
E back to A in (T2"Tj) and using the expression for Co given in (|55p . we get 

A' a 3 \A J (2vr) 4 + P(k) + 0§) 2 J' 1 ] 

Now, using (J[(i 4 /c]£;^A • • • ) = (f[d 4 k]h ■ • • ) — (J[^ 4 ^]s ' ' ' ) an d substituting 
a = fr 1 ^, we can simplify this equation to get 



a 7 " 1 ( + 2 iw (MM) 5 ^ 



wherd^l 



1 /• [d 4 £;] A fc| + l 2 (k) 



A 7 (2tt) 4 (Jfct + Z 2 (fc) + $) 2 ' 



(36) 



12 Note that in this convention, the RG flow is from high to low energies. This is opposite 
to the convention generally used in high energy physics. 

13 In the broken phase, £ = by the gap equation (|20|) . In the unbroken phase, £ is a 
non-zero constant, independent of E. For this reason, it turns out that the RG equation 
obtained in (|38[) . for the dimensionless coupling Aj? defined in (|37p . does not depend on 
£. Consequently the RG equation in the unbroken phase can be obtained from (|3"5)) by 
specializing to 4>r = 0. 
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So, for the dimensionless renormalized coupling, X R = A* 1 A / , we get 



1 1 [ tl0 f [d%E kj + l 2 (k) 



\ R e^ +2 J {ki+i\k)+tl) 2 \ (37) 

This leads to the RG equation 

/ 0- + K) n \9or\ Jo 



l-s 2 (es 2 + /4) 2 

where e = ±1 is the relative sign of go and gi, (pR = (po/ E is the dimensionless 
renormalized coupling corresponding to the vev (with the RG equation cpR = 
<t>n) and 



Note that the form in which the right-hand side of (1381) has been written is 
inappropriate for the special case goR = 0. In this case one must use the 
alternative, but entirely equivalent, form: 

A * = - (~ ~ A * + JriM -irhs I" 1 ds . g2 =,(40) 

/ (1 + 0I) 2 vr 3 |^| 3 7 ^/l _ s 2 ( eyU - 2 s 2 + I) 2 

2 

where h\ = hofi R . It is easy to see that hi — > 1 as /jLr — > oo. 
4.2 The renormalized action 

In terms of the dimensionless renormalized couplings, the low energy action 
( 127)) can be written as 

S = /^^)l(2 7 r) 4 A^o7°-^).7 + 0«V 4 (^-^) 



(2tt) 4 7 (2tt) 4 

7=<r*(* - ?)}**(?) - A 5 " 1 / ^ + Ci**g - C m k^ \a R (k)f 



Ai , — 

— 2— v N(pR<j R (k = 0) + classical part, (41) 
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where l R (k) = k(go R k 2 + gm), k~o = ko/A, k = AA~z. Moreover, the renor- 
malized fields are related to the bare fields by = b^ +l ^i, a R = b^a. 
The coefficients Ci^r are related to and have been defined in ( 160|) . 

This action seems to depend separately on the two couplings g 0R , gm, 
but actually the physics described by it depends only on the combination 
^Ri ( |33|) . To see this more explicitly, let us make the change of variables 
k -> k/\g 1R \, -> \giR\^m, vr -» \gm\ 3 ^R, A^ Ai = \ R /\giR\ 3 . 
After this change of variables, the action takes the form 

S = [ [ ^^R(k){(2^A(hl -lR{k).^ + ^R)s\k- q ) 



(2tt) 4 J (2vr) 4 



A* > — 

-2 — VN(j) R a R (k = 0) + classical part, (42) 
Ai 

where = k(ek 2 / n R + 1) and Ci |2 are given by This form of the 

action makes it explicit that physics depends only on the combination fi R 
since any separate dependence on g 0R , gi R has now disappeared. 

The form (l4"2"j) of the low energy action is not suitable for small values of 
fj,R (equivalently for small values of g± R or large values of goR). In this case, a 
more suitable change of variables in the action (|4T!) is k —* k\go R \~3 , ^i R — > 
\goR\^iR, ctr -> \goR.W R , \r -> A 3 = \ R /\g R.\. After this change of vari- 
ables, the action takes the form 

S = /^^)|(2^) 4 A^o7 O -4(^)-7 + 0hV 4 (^-?) 



(2tt) 4 J (2vr) 4 

]=a R {k - q)}* lR (q) - A^ 1 J^(j- + CjAg - M*)| 5 



Af ^ 

—2—VN(j) R a R (k = 0) + classical part, (43) 
A3 

— » — » — > 2 

where ^(fc) = /c(e/c 2 + fi R ) and C| = Cifi R 2 , C 2 = C^/^ ■ It can be shown 
that C[ 2 have a finite limit as fi R — > 0; see equations (l5Sl) - (f531) . This form of 
the low energy action is now suitable for small values of fi R . 

We have thus found two equally valid descriptions of the physics of the 
4-fermi theory. One is that given by the action in 041|) . which is valid for all 
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the values of the renormalized coupling (j, R . However, in this form the action 
depends on two couplings, g 0R and g 1R . In the form (l4"2j) and (T4"3"j) . the low 
energy action depends only on the combination fi R of these, but two different 
descriptions are needed to cover the entire range of possible values of /j, r . 



4.3 Fixed points 

As we have argued above, the relevant renormalized coupling constants in the 
low energy theory are A x = X R /\gi R \ 3 and A 3 = X R /\g R\, with A 3 = /i|Ai- 
The RG equations for these can be obtained from equations (j35|) and fl4"Ul) 
using ([32]). We get 



(1-0|) A 



2 r h o g 2 



Al = - 2Al + fiT#^/ ds > ■• (45) 



'l - s 2 (efi R 2 s 2 + l) 2 
Together with these, we also have the RG equation for /j, R , namely 

Pr = PR, (Pr 1 )' = Pr 1 - (46) 

The second of these is appropriate for large fi R . Equations pi|) - (Jl"5j) consti- 
tute the set that describes the RG flows in this mode0- We emphasize that 
the explicit dependence on z has dropped out of these equations. This is nice 
since one expects that specific values of z should characterize only the end 
points of an RG trajectory, not the trajectory itself. 

Now, let us first consider the case of small fi R . In this case, the appropriate 
equation is (jUJ). We see that there is a possible fixed point at A3 = 0. For 
this to be a fixed point, we must also have fi R = and (f) R = 0. This is what 
we have been describing as the z = 3 Lifshitz-like fixed point. 

The case /ir — ► 00 is more interesting. In this case we must use (|45|) . 
which in the limit approximates to the equation 



14 Note that 4>r is not an independent variable since it is determined in terms of Ai } 3 
and fi R by the gap equation in the broken phase, while in the unbroken phase it vanishes. 
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For a fixed point we must have (J)r = and one of the two possibilities: 
Ai = 0, 8tt 2 . The first of these is the free field (Gaussian) fixed point and 
the second is a new Lorentz invariance violating fixed point. 

Figure [2] shows a plot of the RG flows near the three fixed points we 
have found. The data for this figure have been obtained using the exact RG 
equations (|44j) and (145|) . Note that we have used e = +1 in these calculations. 
In the broken phase, on the critical line <pR = 0, the RG flow from UV ends 



z = 3 
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Figure 2: In the figure on the left are plotted RG flows in the (A 3 , /ir) plane 
and in the one on the right are flows in the (Ai, I/z^r) plane. fiR increases 
from left to right in both figures. The chiral symmetry broken (unbroken) 
phase is indicated by the letter B (U). The dotted line is the critical curve 
on which 0o vanishes. 



on the Lorentz invariance violating fixed point in the IR. This can be seen 
directly from (1571) . Making the change of variables /c 4 — > Ek^, k — > Ek/\gi\ 
in this equation, we get 

i r jd^u kl + il{k)-<p\ \ 9l \\ 

Ai J (2tt)< (kl + l R (k) + & E*^ 1 1 

where, as before, ln{k) = k(ek 2 /n R + 1) and j[d 4 k]i ■ ■ ■ = J d A k 6(1 — (k\ + 
l%(k))) ■ ■ ■ . Now, in the broken phase, the gap equation implies £ = 0. So, 
for /i^ 1 = 4>r = 0, the right-hand side of the above equation evaluates to 
1/87T 2 . What happens for /j, R = 0, but (J)r ^ 0? In this case, for small 
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values of 4>r, the coupling increases as Ai ~ 87r 2 /(l — 0|,ln0^ 2 ). So, these 
trajectories diverge to larger values of the coupling, doing so faster for larger 
values of (f)R. Figure 2 confirms this for trajectories in the broken phase. 
Note that there is no fixed point for A x — > oo since the beta function of A^ 1 
does not vanish at X^ 1 = 0. 

In the unbroken phase, |(?i| 3 £ is a non-zero constant, independent of the 
flow parameter E, while E — > in the IR. Thus, in the unbroken phase the 
RG flow will terminate at Ai = in the IR. 

What does the theory look like at these two fixed points? Consider first 
the nontrivial fixed point at Ai = 8ir 2 . The fixed point action can be obtained 
from (I42p by setting 0^ = and taking the limit fiR — > oo. For <pR = 0, the 
fermions become massless. In the a kinetic terms, the coefficients Ci,2 grow 
logarithmically in the limit hr — > oo, as shown in ()65p . The implication is 
that as we approach the fixed point Ai = 8n 2 , the kinetic terms for a grow and 
eventually dominate the mass term. This can be seen more directly by the 
rescaling a — > a / \fC\. In the limit llr — > oo, the mass term and the Yukawa 
interaction disappear, leaving behind a massless scalar decoupled from the 
fermions. So the theory at this fixed point has free massless fermions and 
a free massless scalar, the maximum attainable velocity of the latter being 
different from the former, unless Jx 3> 1, in which case Lorentz invariance is 
restored near the fixed point. 

Note that our analysis implies the existence of a fixed point in the usual 
relativistic NJL model as well. This can be established as follows. The 
RG equation for the 4-fermi coupling in the NJL model can be obtained by 
setting p, ^> 1 and z — 1 in (j4"Ul) . Since fiR ^> ft, this also implies [Ir 1. 
For large lir and z — 1, the equation for Ai = A#/|gi#| 3 is precisely ff4T|) . 
Moreover, from the low energy action ( 1421) . we see that in this parameter 
regime the Lorentz violating piece in the fermion kinetic term vanishes and 
the coefficients Ci t2 work out to be those appropriate for a relativistic NJL 
model with a cut-off A, as we have argued below equation (|65|) . 

The other fixed point, that at Ai = 0, is described by just a free massless 
fermion. This is because near this fixed point the a mass goes to infinity as 
Li R because of the manner in which Ai approaches the fixed point, which is 
described by equation ( )48l) . Therefore, this time the rescaling a — > ojyC\ 
leaves the mass term as dominant, with the mass going to infinity. Hence a 
decouples at the fixed point, leaving behind free massless relativistic fermions. 
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This is the theory that one gets from the original four-fermi model at the 
trivial (Gaussian) fixed point. 

5 Concluding remarks 

In this paper we have analysed RG flows in a z = 3 Lifshitz-like four fermi 
model, which is ultraviolet complete in 3 + 1 dimensions. The model flows 
in the infrared to a theory in which Lorentz invariance is violated at 0(1) 
level, which cannot be tuned away by adjusting a parameter. The origin 
of these violations can be traced to fermions with energies higher than the 
Lorentz violating energy scale, propagating in loops and contributing to the 
induced kinetic terms for the composite boson, in the chiral symmetry broken 
vacuum. However, if one works with a finite cut-off, which is taken to be much 
smaller than the Lorentz invariance violating scale, then the model flows in 
the infrared to an approximately Lorentz invariant theory even in the bosonic 
sector, which is similar to the low energy limit of the usual Nambu— Jona- 
Lasinio model in the broken phase. A physical way of interpreting the cut-off 
could be as supersymmetry breaking scale in a supersymmetric version of this 
model. In this case, the offending contributions of fermions in loops would be 
cancelled by their supersymmetric partners and the Lorentz violations would 
be controlled by the ratio of the Lorentz violating scale to the supersymmetry 
breaking scale, which can, in principle, be made small. Possible applications 
of the present model to the Higgs sector of the Standard Model would then 
put constraints on these two scales for consistency with data. 

A remarkable feature of the general RG equations we have obtained, ( 1301) 
and (134)) . is that they describe flow in two scaling parameters, namely a and 
b. More generally, in a theory with anisotropy in n different directions, the 
RG equations will describe flow in n parameters. The parameters presumably 
get related near a fixed point, as in the present example in which we found 
that a = near the fixed point labeled by the exponent z. Away from 
the fixed points, however, a more general flow in multiple scaling parameters 
would seem to be more appropriate. It would be interesting to explore such 
more general flows. 
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A Coefficients appearing in the vacuum so- 
lution 

We give here the coefficients Jo, h, I2 and J3 that appeared in equation (fl8|) : 

h ^ dk TT^ = W^\ 1+ V3 a 2/3 + V31n 2/3 i 

Jo ± + 9o K 12g < V l + ^o ~ 4 V 1 + 9o + V 3# n 



1/3 



J i = / dk ~ — 27^ = — tan 1 9o, 



k 2 J_ 
1 + gfik 6 ~ 3g~ 

d^T, 9i K\9 = ~n — 3 taI1 90 







J 3 = / ^fc - , = f tan^ 1 g + ^ . (49) 

y (i + #) 3 2% 3 V (i + ^o 2 )v 

B Evaluation of the low energy effective ac- 
tion 

Here we give some details of the calculation that lead to the low energy 
action (1271) . Integrating out the high energy modes of the fermions gives the 
effective action 

AS = -i Tr \n(A + B) = -i Tr ln(A) - i Tr ln(l + A~ l B) (50) 

where A(k, q) = (2tt) 4 (7% - 7-f(&) + O ) S 4 (k - q), B(k, q) = ^a(k - q) 
and 'Tr' stands for integration over high momenta (between the lower cut-off 
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E and upper cut-off A, as explained above equation (12 7p ) and trace over all 
indices. AS can be expanded as 

AS = -i Tr ln(A) - i Tr(A _1 5) + ^{A^Bf + higher powers (51) 

Each factor of B comes with a factor of -7= and hence the "higher powers" of 
(A~ l B) are subleading in 1/N. Thus, powers higher than quadratic in u{k) 
in the effective action (!2~Tj) are subleading in 1/N, which we omit. Now, 

Tr \n(A) = iNV j ^fe** tr In (7% - l-l{k) + O ) , (52) 

where 'tr' stands for trace over Dirac indices only. Also, Tr(A~ 1 B) is given 
by 

[d A k} E ->K 1 



Tr(A^B) = Ai(f)Q\fN<j(k = 0) J 



(2tt) 4 k 2 + l\k) + <f) 2 ' 



(53) 



where we have continued the momenta to the Euclidean signature. Thus, the 
coefficient C3 appearing in (j27|) is given by 

r _l o f \d%E^K 1 , . 

3 "A V (2tt)« Jfc2 + l 2 (*)+$" 1 j 



Now, 



(27T) 



xtr 



7% - 7-J(*0 + <h \\ 7°(^o - Po) - 1-Kk - P) + 0o 



Expanding to quadratic order in p , \p\ for small values, the above expression 
gives 

-Tr( A-'B? - I \^k\a(v)\ 2 lav 2 - Cov 2 + 2 I ^ E ~> A k 2 + I 2 (k) - <j> 2 } 

Thus, the action will be of the form given by equation (1271) where Co is given 
by 

r _l 9 f [d%E^ k 2 + l 2 (k)-4> 2 
°"A 2 7 (2*0* (A;i + ^) + ^f (55) 
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Moreover, the coefficients Ci j2 are given by 

OxGto, gu A, E, O ) = / ^|j^{-2A 2 + 4 (/2(A;) + ^ A ' 

-16^(/2 (fc)+ ^ A 4l (56) 



C 2 (^ , <7x; A, E, O ) = | [ ^7 A {A 2 (fc)A 2 - V*) A 3 + A 4 (fc)A 4 |, 
A 2 (fc) = 2((?oA; 2 + 9i)(2g k 2 + ^) + |f7 A; 2 (5f7o^ 2 + 3 gi ), 

A 3 (k) = p 2 (k)(3g k 2 + 9l f + m(g k 2 + 9l )(3g k 2 + <&) 
g 

+ ^go(f>lk 2 (Qg k 2 + A 9l ), 

A 4 (fc) = y 2 / 2 (fc)(3^ 2 + ^i) 2 , (57) 

where A(k 4 , k) = \j(k\ + l 2 {k) + <g). 

For nonzero values of go and gi, one can express the dependence of Ci )2 
on these parameters essentially only through one combination, the scale ^ 
defined in ( 1251) . For example, one can scale out the dependence on g 1 . This 
can be done by the change of the integration variable k — > k/\gi\. We get 

Ci(# ,<7i;A,£,(/>o) = |^r 3 Ci(e/i- 2 ,l;A,E,0o) = \gi\~ 3 C h 
C 2 (<7o,<7i;A,£,^o) = Isil-^e/i-VjA,^) = lyil" 1 ^. (58) 

where e = ±1 is the relative sign of go and g\. In terms of the dimensionless 
renormalized parameters defined in ff29l) and fl33|) . it is easy to show that 

Ci,2 = Cx^e/j,^ 2 , 1; /iif //!, 1, <fi R ) = C 1)2 (e, 1; (59) 

where <f) = <po/A and /i = /i/A are the values of the respective running 
couplings, (fin, /ir, at the UV cut-off. The form in the second equality will 
be useful later when we discuss the E — > limit of these coefficients. Now, 
from the first equality of fl59|) and fl58|) . we get 

Ci = \gin\ 3 Ci(g R,giR] fiR/p>,l,4>R) = \gm\ 3 CiR, 

C 2 = \giR\C 2 (goR, g\R\ vr/P, IAr) = \gm\C2R- (60) 
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Similarly, scaling out the dependence on g , we get 

CiOo,#i;A,£,0 o ) = kr^ita/^A,^^), 
C 2 (2o,<?i;A,£,0 o ) = |^or 1/3 C2(e,//i;A,E,0 o ). (61) 

Moreover, one can easily show that in terms of the renormalized parameters, 

d(e, //§; A, E, O ) = ET^e, 4; 1, = C[/E 2 , 
C 2 (e,^;A,E,cf> ) =£" 2 / 3 C 2 (e,4;^,l,fe) ee C 2 /£ 2 / 3 (62) 
Then, from fl60|) - fl62|) . one gets 

C[ = C 2 = (7 2 /4- (63) 

It is easy to see that C[ 2 are finite for fiR = 0. It follows that the quantities 
on right-hand side of equations in (|6"3"1) have a finite limit as fiR — > 0. 

C Limiting behaviour of Ci 2 

Here we will discuss the limiting behaviour of the coefficients Ci j2 r, given by 
( |56|) . (157|) and (|59|) . for large values of /xr, which is the IR regime of small E. 
We first note that for any positive real number 77 satisfying fl~ 1 > 77" 1 > /x^ 1 , 
we have {j[d%^^ ■ ■ ■ ) = (J[d 4 k] v -^ ■ ■ ■ ) - U[d%^^ •••)• 
It follows from this identity and the second equality of f[5Uj) that 

£1,2 = C 1)2 (e, 1; 1/77, ftp) + C lj2 (e, 1; 1/r,, l/fjL R , $/p) (64) 

We are interested in the limit of large hr for = 0. The fiR dependence of 
C12 comes only from the second term above. In the limit 1 ^> rf x > /i^ 1 , 
the calculation of this term greatly simplifies since the "higher derivative" 
terms can be neglected. That is, in this regime of parameter values, /; C 1 
throughout the integration range, so we may neglect k 2 compared to 1. Then, 
from fl56l) and fl57|) . we find that the leading contribution goes as 

(7 1|2 » C li2 (e, 1; 1/fl, 1/rj, 0) + -Lin . (65) 

We see that Ci j2 grow logarithmically with /i R . If fx ^> 1, then we may choose 
r\ — p. In this case, C\ ~ 6*2 « g^hi (hr/Ji). These are just the values for 
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these coefficients in the relativistic NJL case. However, in general, C\ is 
quite different from (7 2 and both are quite different from the corresponding 
quantities in the relativistic NJL case. Therefore, in general the a kinetic 
terms do not have Lorentz symmetry, as discussed in section 3 below equation 
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